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Abstract: We revisit the impact of the jet algorithm on predictions of energy flow into 
gaps between hard jets, defined using the kt clustering procedure. The resulting predic- 
Qh! tion has two distinct components: a primary emission piece that is related to independent 

emission of soft gluons by the hard jets and a correlated emission (non- global) piece known 
only in the large N c limit. We analytically compute the dependence of the primary emis- 
sion term on the jet algorithm, which gives significantly more insight than our previous 
numerical study of the same. We also point out that the non-global component of the an- 



swer is reduced even more significantly by the clustering than suggested previously in the 
literature. Lastly we provide improved predictions for the latest ZEUS photoproduction 
data, assessing the impact of our latest findings. 
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1. Introduction 

One of the most commonly studied QCD observables is the flow of transverse energy (Et) 
into gaps between jets in various QCD hard processes. Since the Et flow away from jets 
is infrared and collinear safe it is possible to make perturbative predictions for the same, 
which can be compared to experimental data for a given hard process. However since one 
is typically examining configurations where Et is small compared to the hard scale Q of the 
process (e.g. jet transverse momenta in hadronic collisions) the perturbative predictions 
involve large logarithms in the ratio Q/Et. Resummation of logarithmically enhanced terms 
of the form a™ ln n (Q/Et) has proved a challenge that is still to be fully met - complete 
calculations are available only in the large N c limit [jD, [2], 3 . Studies of the E t flow have in 
fact directly led to developments in the theoretical understanding of QCD radiation and 
this process is still ongoing j|. 

Another feature of the energy flow away from jets is its sensitivity to non-perturbative 
effects. Thus one may expect significant 1/Q power corrections to energy flow distribu- 
tions of a similar origin to those extensively studied for various jet-shape observables Q. 
Moreover the E t flow in hadronic collisions is a standard observable used to develop an un- 
derstanding of the underlying event and to assess its role after accounting for perturbatively 
calculable QCD radiation || 0]. 

Given that Et flow studies potentially offer so much valuable information on QCD over 
disparate scales, involving perturbative parameters such as the strong coupling a s , QCD 
evolution, coherence properties of QCD radiation and non-perturbative effects, it is not 
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surprising that they have been the subject of substantial theoretical effort over the past 
few years. 

In this paper we wish to focus on the aspect of resummed predictions for the Et flow 
into gaps between jets. Perhaps the most significant problem involved in making such 
predictions is the non-global nature of the observable [Q, Qj. More precisely in order to 
resum the leading single logarithms involved, one has to address not just a veto on soft 
emissions coupled to the underlying primary hard parton antennae (known as the primary 
emission term), but additionally correlated emission or non- global contributions, where a 
clump of energy-ordered soft gluons coherently emit a still softer gluon into the gap region 
O. For this latter contribution the highly non-trivial colour structure of the multi-gluon 
"hedgehog" configuration has proved at present too significant an obstacle to overcome. 
One thus has to resort to the large iV c limit to provide merely a leading-log estimate for 
the away-from-jet Et flow. This situation can be contrasted with the case of event-shapes 
and Drell-Yan qx resummations which have been pushed to next-to-leading and next- 
to-next-to-leading logarithmic accuracy respectively. The impact of finite N c corrections 
in non-global observables is thus a factor in the theoretical uncertainty involved in the 
corresponding resummed predictions. 

Given that the non-global component has a substantial quantitative impact over a 
significant range of Et values for a given hard scale Q and that it is computable only 
in the large N c approximation, it is clearly desirable to reduce the sensitivity of a given 
observable to non-global logarithms. An important observation in this regard was made in 
Ref. ||: if one employs the kt clustering algorithm 0, 10 1 to define the final state such that 



the energy flow into a gap between jets is due to soft /ct-clustered mini-jets (rather than 
individual hadrons), the non-global logarithms are significantly reduced in magnitude 1 . 
This observation was exploited to study the case of Et flow in dijet photoproduction where 
a result was provided for the primary emission component of the Et distribution and the 
reduced non-global component was modeled |T^| . 

However it has subsequently been found that kt clustering also has a non-trivial impact 
on the primary emission component of the result [13]. This was not taken into account in 



Refs. [||, 12] and also affects the ability to make resummed predictions for a host of other 
jet observables such as azimuthal correlation between jets ^(pjj- I n f ac t the findings of 
Ref. [13] are not just specific to the kt algorithm but would also crop up in the case of jet 
observables defined using iterative cone algorithms. 

In the present paper we wish to shed more light on the resummation of the primary or 
independent emission component of the result and its dependence on the clustering algo- 
rithm. While the leading O (a 2 s \v? (Q / Et)) clustering dependent behaviour was computed 
analytically in Ref. [13], the full resummed result for the primary emission component was 



computed only numerically in the case of a single hard emitting dipole (e + e~ — > 2 jets 
or DIS 1 + 1 jets). Here while sticking to a single hard dipole we shed more light on the 
structure of the primary emission term and analytically compute it to an accuracy that is 
sufficient for a wide range of phenomenological applications. 



1 For recent progress on aspects of the kt algorithm itself see Ref. [ill . 
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The analytical insight and calculations we provide here will also make the generalisation 
of the /^-clustered primary emission result to the case of several hard emitters (dijets 
produced in photoproduction or hadron-hadron processes), involving a non-trivial colour 
flow, relatively straightforward. 

The above resummation of the primary component of the answer assumes greater 
significance when we discuss our second observation: once an error is corrected in the 
numerical code used for the purposes of Refs. & |l2| the non-global component of the 
result is reduced even more compared to the earlier estimate. With a very small non- 
global component (which can be numerically computed in the large N c limit) and a primary 
emission component that correctly treats the dependence on the jet algorithm, one is better 
placed to make more accurate resummed predictions than has been the case till now. This 
is true not just for the E t flow but also as we mentioned for a variety of jet observables 
for which there are either no resummed predictions as yet, or only those employing jet 
algorithms not directly used in experimental studies jL4|. 

This paper is organised as follows. In the following section we define the observable in 
question and revisit the issue of the dependence of the primary and non-global pieces on the 
jet clustering algorithm. Following this we demonstrate how the primary or independent 
emission piece can be computed at all orders in a s , accounting to sufficient accuracy for 
the effects of the clustering algorithm. We explicitly describe the case of three and four- 
gluon contributions to demonstrate the steps leading to our all-order results. Following this 
we re-examine the non-global component of the answer and find that this is significantly 
smaller than earlier calculations of the same ||] . We put our findings together to examine 
their impact on photoproduction data from the ZEUS collaboration [[0|] and lastly point 
to the conclusions one can draw and future extensions of our work. 

2. Resummation of the primary emissions 

Let us consider for simplicity the process e + e~ — > 2 jets. The calculations for processes 
involving a larger number of jets and more complex jet topologies can be done along similar 
lines. 

We wish to examine the Et flow in a region Q which we choose as a slice in rapidity 2 
of width Ar] which we can centre on 77 = 0. We then define the gap transverse energy as: 



where the index i refers to soft jets obtained after kt clustering of the final state. We shall 
concentrate on the integrated E t cross-section which is defined as: 






(2.2) 



with a the total cross-section for e + e — > hadrons, with center-of-mass energy Q. 



2 Since we are here dealing with back-to-back jets we can define the rapidity with respect to the jet axis 
or equivalently, for our purposes, the thrust axis. 
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The single-logarithmic result for the above, without kt clustering (where the sum over 
i in Eq. (2.1) refers to hadrons in the gap rather than jet clusters), was computed in Ref. Q 



and can be expressed as: 

E(Q, Qn) = E P (t) S(t) , t = - -± a s (k t ) . (2.3) 

2?r J Qn h 

The above result contains a primary emission or "Sudakov" term 3 £p(t) and a non-global 
term S(t). 

The primary emission piece is built up by considering only emissions attached to the 
primary hard partons namely those emitted from the hard initiating qq dipole in our 
example, while the non-global term arises from coherent soft emission from a complex 
ensemble of soft emitters alongside the hard initiating dipole. More precisely we have: 

E P (t) = e - ACFtAr > , (2.4) 

which is the result of resumming uncancelled fcj-ordered virtual-emission contributions, in 
the gap region. The non-global component, as we stated before, is computed numerically 
in the large N c limit. 

Next we turn to the fet-clustered case. The result stated in Ref. || assumes that the 
primary or Sudakov piece is left unchanged by clustering since it appears to be the expo- 
nentiation of a single gluon emitted inside the gap. The non-global piece is recomputed 



numerically implementing clustering ||. As already shown in Ref. [ 13 1 however, the as- 
sumption regarding the primary emission piece being unaffected is in fact untrue and this 
too needs to be recomputed in the presence of clustering. The corrections to the primary 
emission term first appear while considering two gluons emitted by the hard qq dipole 
and persist at all orders. Below we provide a reminder of the two-gluon case discussed 



in Ref. 13] and subsequently consider explicitly the three and four-gluon emission cases 



before writing down the result to all orders as a function of the radius parameter R. 
2.1 Two-gluon emission 

In order to examine the role of the kt algorithm we point out that in our case (fc^-ordered 
soft limit) one can start the clustering procedure with the lowest transverse-energy parton 
or equivalently the softest parton. One examines the "distances" of this particle, i, from its 
neighbours, defined by dij = E 2 ^ ^(Ar^) 2 + (A^) 2 ^, where E t) i is the transverse energy 
of the softest parton. If the smallest of these distances is less than E^R 2 , particle i is 
recombined or clustered into its nearest neighbour and the algorithm is iterated. On the 
other hand if all dij are greater than E^R 2 , i is counted as a jet and removed from the 
process of further clustering. The process continues until the entire final-state is made 
up of jets. Also in the limit of strong energy-ordering, which is sufficient to obtain the 
leading-logarithms we are concerned with here, the recombination of a softer particle with 
a harder one gives a jet that is aligned along the harder particle. 



3 We use the term "Sudakov" in a loose sense since the primary emission result leads to an exponential 
that is analogous to a Sudakov form-factor. 
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The dependence of the primary emission term on the jet clustering algorithm starts 
naturally enough from the two-gluon level. While the Sudakov result exp (— ACptArj) comes 
about due to assuming real-virtual cancellations such that one is left with only virtual 
emissions with kt > Qn m the gap region (for the integrated distribution), Jq clustering 
spoils this assumed cancellation. 

Specifically let us take two real gluons k\ and &2 that are ordered in energy {uj\ S> w?). 
We consider as in Ref. [13] the region where the softer gluon kz is in the gap whilst the 



harder k% is outside. Additionally we take the case that the gluons are clustered by the 
jet algorithm which happens when (At/) 2 + (Acj)) 2 < R 2 with Arj = r\% — r\\ and similarly 
for A(p, which condition we shall denote with the symbol #21- Since k-z is clustered to 
k\ it gets pulled outside the gap, the recombined jet being essentially along k\. Thus in 
this region the double real-emission term does not contribute to the gap energy differential 
distribution da/dE t . Now let k\ be a virtual gluon. In this case it cannot cluster k2 out 
of the gap and we do get a contribution to the gap energy differential distribution. Thus a 
real-virtual cancellation which occurs for the unclustered case fails here and the mismatch 
for the integrated quantity E(i), amounts to: 



cr thrill drjt— I ^fa-HQffi.g, (2.5) 



where we reported above the result computed for R < Ary, in Ref. [|T^]. Here we intro- 
duced the primary emission term Cn that corrects the Sudakov result at 0(a™) due to the 
clustering requirement. 

The fact that the result scales as the third power of the jet radius parameter is interest- 
ing in that by choosing a sufficiently small value of R one may hope to virtually eliminate 
this piece and thus the identification of the primary result with the Sudakov exponent 
would be at least numerically accurate. However the non-global term would then be signif- 
icant which defeats the main use of clustering. If one chooses to minimise the non-global 
component by choosing e.g. R = 1, then one must examine the primary emission terms in 
higher orders in order to estimate their role. To this end we start by looking at the three 
and four-gluon cases below. 

2.2 Three-gluon emission 

Consider the emission of three energy-ordered gluons fci, k<± and k% with < W2 <C u)\, 
off the primary qq dipole, and employing the inclusive kt clustering algorithm 10 as 
explained previously. 

We consider all the various cases that arise when the gluons (which could be real or 
virtual) are in the gap region or outside. We also consider all the configurations in which 
the gluons are affected by the clustering algorithm. We then look for all contributions where 
a real-virtual mismatch appears due to clustering, that is not included in the exponential 
Sudakov term. The Sudakov itself is built up by integrating just virtual gluons in the gap, 
above the scale Qq. The corrections to this are summarised in table |]. 

In order to obtain the various entries of the table one just looks at the angular con- 
figuration in question, draws all possible real and virtual contributions and looks for a 
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Table 1: Contributions of different configurations of particles to Sp(t) at O(al). We define 
% =8(R 2 - (r/i - 7]j) 2 - (fa - fa) 2 ), e.g. 9 13 = 1 means (rji - r/ 3 ) 2 + (fa - fa) 2 < R 2 . We also 
define W = (— 4Cpi) 3 /3!, so the entries "W" indicate a miscancellation which leads to a single-log 
correction to the Sudakov result, while the entries "0" indicate a complete real-virtual cancellation. 
We have discarded the case where all particles are in the gap since such configurations are already 
included in the exponential Sudakov result. 



mismatch between them generated by the action of clustering. We translate table |l| to: 



x{ / d m^p- [ dm^r- / d7] 3 e 32 9 31 + 

l Jk^n Ait Jk 2 ^Q Z7r Jk :i en 

+ [ dffjp- I d V2 ^ f d V3 [#31 + (1 - 03l)(l - M#2l] + 

Jk^n Ait J k2 £Q lit Jk s en 

+ f drji — — f drj2 — [ drj 3 832 \ , 
Jkren Ait J k2 £Q lit J ka€ Q J 



where we used the freedom to set 03 = 0. We identify three equal contributions consisting 
of the integrals in which there is only one theta function constraining only two particles: 
the last integral over #32, the integral over #31 and that over #21 in the third line. The 
set of configurations #32, #31 and #21 is just the set of constraints on all possible pairs of 
gluons, and in fact we can generalise the factor 3 to the case of any number n of gluons by 
n(n — l)/2, which will enable us to resum i? 3 terms. We shall return to this observation 
later. The integrals of the above type reduce essentially to the clustered two-gluon case as 
calculated in Eq. fl2.5|), and the integral over the third "unconstrained" gluon is just Ar]. 
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Explicitly we write Eq. as: 



dvi^r- / d7] 2 — / dr/3 6*32 6/31 + 
+ / d Vl ^ [ [ dTfe [0310KJ -031 -032] 021 + 

+3x/ c^/ d^/" d^flsaV (2.7) 



Computing the various integrals above (for simplicity we take -R < Ar]/2, which is sufficient 
for our phenomenological purposes) one obtains: 



ir 32 \ R 5 i? 5 /vr 32 \ R 5 32 R 5 2 A . . 



with / ~ 0.2807 and we have written the results in the same order as the five integrals 
that arise from the various terms in Eq. fl2.7|). Hence: 



Cf = -^CptY l 3 x —At? R* + h R b } , (2-9) 



where f% — —0.04360. We note the appearance of an R 5 term which, as we shall presently 
see, persists at higher orders. This term is related to a clustering constraint on three gluons 
at a time via the product of step functions 632 #21 (#31 — 1) with ii2, ^3 £ ^ and k\ ^ fl. 

Next we look at the emission of four soft, real or virtual energy-ordered gluons. This 
will help us move to a generalisation with any number of gluons. 



2.3 Four-gluon case and beyond 

Now we take the case of four-gluon emission and identify the patterns that appear at all 
orders. A table corresponding to table p] is too lengthy to present here. The result can 
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however be expressed in an equation similar to that for the three-gluon case. We have: 



x< / / / 43 + 

'lin J2m .73 out J4in 



+ / / / / [#42 + ^32(1 M] + 

J 1 in J 2 out J 3 in J 4 in 

+ / / / / {#41 + 0-41 [031 0-43 + #43 #21 0-42 + #21 0-42 0-43 0-31^-32]} 
J lout J 2 in J 3 in J 4 in 



+ 



+ / / / / 042 043 + 
'lin J 2 out J 3 out J4in 



+ 



+ / / / / 043 [041 + _41 0-42 02l] + 
J 1 out J 2 in J 3 out 1/ 4 in 

+ / / / / 041 042 + 041 0-42 0-43 032 + 0-41 0-43 031 [042 + 0-42 032] 

J 1 out J 2 out J 3 in J 4 in 

+ f I f I 041 042 04 3 ), (2.10) 
J lout J 2 out J 3 out J 4 in J 

where 0_jj = 1 — 0jj and "in" or "out" pertains to whether the gluon is inside the gap 
region or out. For brevity we did not write the differential phase-space factor for each 
gluon which is as always drj d(j>/(2ir). We identify six R 3 terms exactly of the same kind as 
computed before and similarly four R 5 terms. Explicitly we have: 

CI = ^{-±C F tf x 



X<|6 X / / / / 043 + 

'lin J2in J 3 out J4in 



+4 X 



1 in J 2 out J 3 out J 4 in 



042 043 + / / / / 032 [043 042 ~ 043 ~ 042] ) + 

J l in J 2 out J 3 in J 4 in / 

+3X ! ft I 21 43 [1 " 041 " 042 + 041 042] + 
J 1 out J 2 in J 3 out J 4 in 

+ / / / / 021 [042 043 — 042 ~ 043 _ 041 0-42 0-43] [031 032 ~ 031 _ 032] + 

J 1 out J 2 in J 3 in J 4 in 

+ / / / / 032 031 [04l(l " 04 3 )(042 " 2) - 43 ] + 
J 1 out J 2 out J 3 in J 4 in 

+ l l l l #41 042 043 )• (2.H) 
J 1 out J 2 out J 3 out J 4 in J 

We discuss below each set of integrals, generalise the result to the case of n emitted gluons 
and then resum all orders. 

• The integral: 

h-4C F t) 4 6x III I 043- (2.12) 

^- Jlin J2in J3out J4in 
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The integrals over particles 1 and 2 give (A77) 2 . The remaining integrals reduce to the result 
computed for the two-gluon case, i.e. the R? term, multiplied by a factor of 6 accounting 
for the number of pairs of gluons n(n — l)/2, for n = 4. Explicitly we have for this term: 

^-4C F t)^Ar,^R* . (2.13) 

For n emitted gluons the i? 3 term, which is always related to the clustering of two gluons, 
is given by: 

i^l-^V^- ( 2 - 14 ) 

TV. 2 D7T 

Hence to all orders one can sum the above to obtain: 

e -4C F tA v (-4CFt) 2 2_ r3 (215) 



• The integrals: 

l(-4C F t) 4 4x ( [ III #42#43 + 
4! V Jlin J2out J3out Jim. 

+ 11 II °32 [#43 #42 - #43 ~ #4 2 ] ) ■ (2.16) 
J 1 in J 2 out J 3 in J 4 in / 

The integral over particle 1 gives Arj, while the rest of the integrals reduce to the ones 
calculated earlier which gave the R 5 result, accompanied with a factor of 4 standing for the 
number of triplet combinations formed by four gluons. For n emitted gluons this factor is 
n(n — l)(n — 2)/3!. Explicitly we have for this case: 



^(-4C F ^^^A V ^f 2 R 5 . (2.17) 



At the n th order we obtain: 



±(-4C F tA V r n{n ^ 2) ArT 3 /2fl 5 , n>3. (2.18) 



Summing all orders we get: 



-iC F tA v (-4C F t) 3 hR5 (21Q) 



• The integral: 



i(-4C F t) 4 3x I I I I #2i.#.i: i . (2.20) 



1 out J 2 in J 3 out J 4 in 



This integral can be factored into two separate integrals involving the constraint on k± and 
&2 and over £3 and k 4 respectively. Each of these reduces to the R 3 result obtained in the 
two-gluon case. Thus we get: 

l -(-4C F t) 4 3x (—) R 6 . (2.21) 



4! v ' \3tt 
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At n th order this becomes: 



ln(n-l)(n-2)(n-3), ^ m \«a -4/^ 



2 



, 8 (-4C F tA V ) n Ar/- 4 ^— j n>4, (2.22) 

which can be resummed to: 

e -4C Ff A, (-^) 4 (|^) 2 ^ 6 - ( 2 - 2 3) 

The factor 3 (and generally n(n — — 2)(n — 3)/8) is the number of configurations formed 
by four (and generally n) gluons such that we have two pairs of gluons each is formed by 
an out-of-gap gluon connected to a softer in-gap one. 

• The remaining integrals 

These integrals give at most an 0(R 7 ) term because they constrain all the four gluons at 
once. In fact for gap sizes A77 > 3R, these integrals go purely as R 7 with no dependence 
on A77. Since here however we wish to use the condition Arj > 2R, which allows us to 
make use of the whole range of HERA data, these integrals do not depend purely on R 
but are a function of R and Ar] which have an upper bound of order R 7 . This can be seen 
by noting that there are three azimuthal integrations that each produce a function which 
has a maximum value proportional to R, so the result of integrating over all azimuthal 
variables is a factor that is bounded from above by i? 3 . Similarly there are four rapidity 
integrations with a clustering constraint on all four gluons implying that they can produce 
an R 4 term at most. In general the result at n th order of constraining n gluons at once, is 
bounded from above by a factor of order R 2n ~ l . 

We can write the result for all these as (— 4Ci?t) 4 /4! y(R, Ar/), and resum such terms 
to all orders (in the same manner as before) to: 

e ^tM yWA|j)) (2 . 24) 

where y(R,Arj) is at most 0(R 7 ). We do not calculate these terms (though it is possible 
to do so) since the accuracy we achieve by retaining the R 3 , R 5 and R 6 terms, we have 
already computed, is sufficient as we shall show. 

The five-gluon case is too lengthy to analyse here. The same patterns as pointed out 
above persist here but new terms that are at most 0(R 9 ) appear when all five gluons are 
constrained. There is also an R 8 term, coming from the combination of R 3 and i? 5 terms 
in the same manner that the R 6 term arose as a combination of two R 3 terms. 



3. All-orders result 

From the above observations we can assemble an all-orders result to R e accuracy, where 
we shall consider R to be at most equal to unity. The final result for primary emissions 
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alone and including the usual Sudakov logarithms (for Arj > 2R) is: 
E P (t) = e - 4C ^x 

1 + (-4C F t) 2 ^R 3 + (-4C^) 3 ^# 5 + {-ICptf-^R 6 + ( ~ 4 ^ F * )4 0(i? 7 ) 

(3.1) 

Formally one may wish to extend this accuracy by computing a few more terms such 
as those integrals that directly give or are bounded by an R 7 behaviour and this is possible 
though cumbersome. It should also be unnecessary from a practical viewpoint, especially 
keeping in mind that R = 0.7 is a preferable value to R = 1, in the important case of 
hadron collisions 4 and the fact that even at R = 1 the R? term significantly dominates the 
result over the range of t values of phenomenological interest, as we shall see below. 

We further note that if one keeps track of all the terms that come about as a combina- 
tion of ill 3 and R 5 terms in all possible ways at all orders, one ends up with the following 



form for Eq. (3.1): 



,(,) _ exp £tf + £Wt h ^ + HgaV)) , (3 . 2) 



the expansion of which agrees with Eq. ( |3.1| ). In the above by 0{R 7 ) we mean terms that, 
while they may depend on Arj, are at most as significant as an R 7 term. We also mention 
that in the formal limit Arj — > oo, there is no dependence of the clustering terms on Ar] 
and they are a pure power series in R. The limit of an infinite gap appears in calculations 
where the region considered includes one of the hard emitting partons. An example of such 
cases (which have a leading double-logarithmic behaviour) is once again the quantity A<j)jj 
between jets in e.g. DIS or hadron collisions. 

Fig. H represents a comparison between the leading R 3 result (i.e. the pure fixed-order 
result of Ref. |l3| combined with the resummed Sudakov exponent), the resummed -R 3 , 
i? 5 and -R 6 result (Eq. ( |3.1| )) and a numerical Monte Carlo estimate with and without 
clustering. The Monte Carlo program in question is essentially that described in Ref. 
with the modification of kt clustering where we computed just emissions off the primary 
dipole "switching off" the non-global correlated emission. 

We note that the resummed analytical form ( |3.1| ) is in excellent agreement with the 
numerical result which contains the full R dependence. We have tested this agreement 
with a range of values of R. We take this agreement as indicating that uncomputed R 7 
and higher terms can safely be ignored even at R = 1 and even more so at fractional values 
of R, e.g. R = 0.7. To provide an idea about the relative role of terms at different powers 
of R in Eq. ( |3,1| ) we note that for R = 1 and t = 0.25 the resummed it! 3 term increases the 
Sudakov result exp (— ACptArj) by 19%, the i? 5 term represents a further increase of 1.5% 
to the result after inclusion of the resummed R 5 term and the R e term has a similar effect 
on the result obtained after including up to R 5 terms. 

Next we comment on the size of the non-global component at different values of R. 



4 This is because the underlying event will contaminate jets less if one chooses a smaller R. 
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Figure 1: Comparison of the analytical results to a numerical Monte Carlo estimate. 
4. Revisiting the non-global contribution 

We have seen above how the primary emission piece is dependent on the jet clustering 
algorithm. It was already noted in Ref. || that the non-global contribution is significantly 
reduced by clustering. Here we wish to point out that after correction of an oversight in 
the code used there, the non-global component is even more significantly reduced than 
previously stated in the literature. Indeed for R = 1 and the illustrative value of t = 0.15, 
which corresponds to gap energy = 1 GeV for a hard scale Q = 100 GeV, the non-global 
logarithms are merely a 5% effect as opposed to the 20% reported previously || and the 
over 65% effect in the unclustered case. 

In Fig. [2] we plot the curves for the primary and full results (in the large N c limit) for 
the integrated quantity E(t) as a function of t defined earlier. We note that for R = 0.5 the 
primary result is essentially identical to the Sudakov result. The non-global contribution 
(which is the ratio of the full and primary curves) is however still quite significant. Ne- 
glecting it leads to an overestimate of 40% for t = 0.15. Increasing the jet radius in a bid 
to lower the non-global component we note that for R = 0.7 the impact of the non-global 
component is now just over 20% while the difference between the full primary result and 
the Sudakov result is small (less than 5%). The situation for R = 1 is a bit different. Here 
it is the non-global logarithms that are only a 5% effect (compared to the 20% claimed 
earlier ||) while the full primary result is bigger than the Sudakov term by around 11%. 

The value R = 1 is in fact the one used in the HERA analyses of gaps-between-jets 
in photoproduction. It is now clear that such analyses will have a very small non-global 
component and a moderate effect on primary emissions due to clustering. In order to 
completely account for the primary emission case for dijet photoproduction one would 



- 12 - 



need to generalise the calculations presented here for a single qq dipole to the case of 
several hard emitting dipoles. An exactly similar calculation would be needed for the case 
of hadron-hadron collisions and this is work in progress. It is straightforward however to 
at least estimate the effect of our findings on the photoproduction case and we deal with 
this in the following section. 



5. Gaps between jets at HERA — the ZEUS analysis 

We can test the perturbative framework presented in this paper with energy flow measure- 
ments in the photoproduction of dijets. These energy flow observables are defined with 
two hard jets in the central detector region separated by a gap in pseudorapidity. A gap 
event is defined when the sum of the hadronic transverse energy in the gap is less than 
a cut-off, and the gap fraction is defined as the ratio of the gap cross-section to the total 



inclusive cross-section. The energy flow observables measured by ZEUS [15] and HI [16] 
use the kt clustering definition of the hadronic final state, and the transverse energy in 
the gap is given by the sum of the mini-jet contributions. In this paper we focus on the 
ZEUS measurements and provide revised theoretical estimates for them. These revisions 
lead to changes that are minor in the context of the overall theoretical uncertainty but 
should become more significant once the matching to fixed higher-orders is carried out and 
an estimate of the next-to-leading logarithms is obtained. The HI data was considered in 
Ref. |TI|, where the theoretical analysis consisted of only the resummed primary emission 
contribution without taking account of the effect of k t clustering. 

The ZEUS data was obtained by colliding 27.5 GeV positrons with 820 GeV protons, 
with a total integrated luminosity of 38.6 ±1.6 pb _1 in the 1996-1997 HERA running 



period. The full details of the ZEUS analysis can be found in Ref. |15], but the cuts 
relevant to the calculations in this paper are: 

0.2 < y < 0.75, 

Q 2 < 1 GeV 2 , 
6, 5 GeV < E T (1,2) , 

|t/(1,2)| <2.4, 
|0.5(r/i + 77 2 )| < 0.75, 

2.5 < Arj < 4, 

where y is the inelasticity, Q 2 is the virtuality of the photon, £^(1,2) are the transverse 
energies of the two hard jets, 7/(1, 2) are the pseudorapidities of the two hardest jets and A77 
is the jet rapidity difference. The further requirement for the gap sample is £t, ga p < Qn = 
0.5, 1, 1.5, 2 GeV, and the clustering parameter R is always taken to be unity. 

The theoretical prediction for the gap fraction is composed of the primary piece, with 
corrections due to clustering, and the non-global piece. We shall now describe each in turn. 

The resummed primary contribution ignoring the clustering corrections, is obtained 



from the factorisation methods of Sterman et al [14] and is described in Ref. 12]. The 



four-jet case of photoproduction requires a matrix formalism, and the exponents of the 
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Sudakov factors in the gap cross-section are anomalous dimension matrices over the basis 
of possible colour flows of the hard sub-process. The emission of soft gluons cause mixing 
of the colour basis. Consideration of the eigenvectors and eigenvalues of the anomalous 
dimension matrices, together with sub-process-dependent hard and soft matrices, allows 
the resummed four-jet primary emission differential cross-section to be written as [12]: 



% = Y.H IL s LI ex P ! [ (xi( v ,n) + x^n)) ° ^a.(/*)j , (5.1) 



where H and 5 denote the hard and soft matrices (expanded over the colour basis), A 
denotes the eigenvalues of the anomalous dimension matrices, 77 = Ary/2 and pt is the hard 



scale of the process. This was computed in Ref. |12] for the case of photoproduction and 
energy flow observables measured by HI. In this paper we have recomputed this differential 
gap cross-section for the observable defined by the ZEUS collaboration. The uncertainty 
in the renormalisation scale is quantified by varying the hard scale in the resummation by 
a factor of 2 (upper bound) and 0.5 (lower bound). 

We now need to account for the effect of clustering on Eq. (^l|). Since we do not 
have as yet the full results for the four-jet case of photoproduction we simply estimate the 
full correction as the square of the correction arising in the two-jet case dealt with here, 
using the appropriate colour factors for each hard sub-process. This was also the method 
used to approximate the non-global contribution for the four-jet case in Ref. [|E|. While 
we emphasise that this is only a rough way of examining the impact of the clustering 
dependent terms computed here, given the size of the effects we are dealing with, it is 
clear that no significant differences ought to emerge if one were to properly compute the 
various dipole sub-processes we need to account for. We also include the revised and 
virtually negligible non-global component in an identical fashion to arrive at the best 
current theoretical estimates. 

The results for the ZEUS gap- fraction with a /ct-defined final state are shown in Figs. || 
and ||. We consider here two different values for the gap energy Qq. For the value of 
Qq = 0.5 GeV one notes that the full prediction accounting approximately for all additional 
sources of single-logarithmic enhancements, is somewhat higher than the pure "Sudakov" 
type prediction. This is due to the extra primary terms we compute here, non-global 
corrections being negligible. For a larger value of = 1.0 GeV the difference between the 
clustered and unclustered primary results is negligible. We also note the large theoretical 
uncertainty on the prediction as represented by the renormalisation scale dependence. 
This is to be expected in light of the fact that the predictions here are not matched to 
fixed-order and account only for the leading logarithms. Improvements along both these 
directions should be possible in the immediate future after which the role of the various 
effects we highlighted here should be revisited. 

6. Conclusions 

In the present paper we have shed further light on resummations of fct-clustered final states. 
We have shown that both the primary and non-global components of the resummed result 
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are affected by clustering and dealt with the resummation of each in turn. For the non- 
global component we find that the results after applying clustering are different from those 
presented earlier ||. The new results we present here indicate an even smaller non-global 
component than previously believed. 

We have also shown how the primary emission clustering effects can be resummed to 
all orders as an expansion in the clustering parameter R and computed a few terms of the 
series. The analytical results we have provided here for a single emitting dipole should be 
generalisable to the case of several hard dipoles (multi-jet processes). This should then 
enable one to write a correct resummed result for primary emissions to a high accuracy 
and deal with the reduced non-global component in the large iV c limit. Such progress is 
relevant not just to energy-flow studies but to any jet observable of a non-global nature, 
requiring resummation. An example is the azimuthal angle A</>,-,- between jets, mentioned 
previously. The work we have carried out should enable next-to-leading log calculations of 
such jet observables to sufficient accuracy to enable phenomenological studies of the same. 

Lastly we have also mentioned the impact of the new findings on the ZEUS gaps- 
between-jets analysis. Since the non-global effects are very small for R = 1 the main 
new effect is the additional clustering dependent primary terms we computed here. Ap- 
proximating the effect of these terms for the case of photoproduction, somewhat changes 
the theoretical predictions but this change is insignificant given the large theoretical un- 
certainty that arises due to missing higher orders and unaccounted for next-to-leading 
logarithms. We consider both these areas as avenues for further work and hope that more 
stringent comparisons can thus be made in the very near future. 
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Figure 2: Comparison of the Sudakov result, the correct primary result and the full result including 
non-global logarithms, for different values of R and with Arj = 1. All quantities are shown in the 
large N c limit for ease of comparison. 
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Figure 3: The gap fraction for the ZEUS analysis with a ^-defined final state (R = 1.0 and 
Qn = 0.5 GeV). The solid line shows the effect of resummed primary emission, the primary emission 
clustering correction factor and the non-global suppression factor. The overall theoretical uncer- 
tainty in all three contributions is shown by the dotted lines. The dashed line indicates the gap 
fraction obtained by only including primary resummed emission without accounting for clustering. 
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Figure 4: The gap fraction for the ZEUS analysis with a kt-defmed final state (R = 1.0 and 
Qn = 1.0 GcV). The solid line shows the effect of resummed primary emission, the primary emission 
clustering correction factor and the non-global suppression factor. The overall theoretical uncer- 
tainty in all three contributions is shown by the dotted lines. The dashed line indicates the gap 
fraction obtained by only including primary resummed emission without accounting for clustering. 
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